A/-A(1232) axial form factors from weak pion production 
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The NA axial form factors are determined from neutrino induced pion production ANL & BNL 
data by using a state of the art theoretical model, which accounts both for background mechanisms 
and deuteron effects. We find violations of the off diagonal Goldberger-Treiman relation at the level 
of 2a which might have an impact in background calculations for T2K and MiniBooNE low energy 
neutrino oscillation precision experiments. 
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The A(1232) resonance is the lightest baryonic excita- 
tion of the nucleon. In addition, it couples very strongly 
to the lightest meson, the pion, and to the photon. As 
a consequence, the A(1232) is of the utmost importance 
in the description of a wide range of hadronic and nu- 
clear phenomenology going from low and intermediate 
energy processes ll|, |2J to the GZK cut-off of the cos- 
mic ray flux 0, |j|. On the other hand, despite its 
large width, it is well separated from other resonances 
what facilitates its experimental investigation. In partic- 
ular, the electromagnetic nucleon to A(1232) excitation 
processes, induced by electrons and photons, have been 
extensively studied at many experimental facilities like 
LEGS, BATES, ELSA, MAMI, and J-LAB. For a recent 
review see Ref. [5|, where also many of the recent the- 
oretical advances in the understanding of the resonance 
have been addressed. 

There has also been a great theoretical interest in the 
axial nucleon A transition form factors. Recently, they 
have been studied using quark models [6|, Light Cone 
QCD Sum Rules 0, Lattice QCD and Chiral Pertur- 
bation Theory (%PT) @, lioj ]. These form factors are of 
topical importance in the background analysis of some of 
the neutrino oscillation experiments (e.g. However, 
their experimental knowledge is less than satisfactory. 
Although the feasibility of their extraction in parity- 
violating electron scattering has been considered 12j , the 
best available information comes from old bubble cham- 



ber neutrino scattering experiments at ANL [13j, |14| and 
BNL 15|, |16|. These experiments measured pion pro- 



duction in deuterium at relatively low energies where the 
dominant contribution is given by the A pole (AP) mech- 
anism: weak excitation of the A(1232) resonance and its 
subsequent decay into Nw. Only very recently, 7r° pro- 
duction cross sections have been measured at low neu- 
trino energies and with good statistics |17j . However, 



the target was mineral oil what implies large and diffi- 
cult to disentangle nuclear effects. Thus, these data are 



less well suited for the extraction of the NA axial form 
factors. 

Besides the original experimental publications, there 
are many studies of the ANL and/or the BNL data in 
the literature [H-|23| with different advantages and short- 
comings. Some of those studies are discussed below. In 
this letter, we analyze the ANL and BNL data incorpo- 
rating the deuteron effects, with a proper consideration of 
statistical and systematical uncertainties and taking ad- 
vantage of several recent developments: improved vector 
form factors and a new model for weak pion production 
off the nucleon that includes background terms. 



A convenient parameterization of the 



ver- 



tex is given in terms of eight q (momentum transfer 
square) dependent form-factors: four vector and four ax- 
ial (C^4 5 6 ) ones. We follow the conventions and no- 
tation of Ref. 



21 



Vector form factors have been deter- 
mined from the analysis of photo and electro-production 
data. Here, we use the parameterization of Lalakulich 
et al. [24], as done in Ref. 21 1. Among the axial form 



factors the most important contribution comes from C5 . 
The form factor Cg , which contribution to the differen- 
tial cross section vanishes for massless leptons, can be 
related to C§ thanks to the partial conservation of the 
axial current (Cg^g 2 ) = C$(q 2 )- 
M the pion and nucleon masses, 



, with 
respectively) 



m x and 
Since 



there are no other theoretical constraints for C^ 45 (g 2 ), 
they have to be fitted to data. Most analysis, includ- 
ing the ANL and BNL ones, adopt Adler's model (25j 
where 1 C^(q 2 ) = and C£ (g 2 ) = -C^(g 2 )/4. For 
C§ several q 2 parameterizations have been used [3, [22| , 
though given the limited range of statistically signifi- 
cant q 2 values accessible in the ANL and BNL data, it 



1 Setting to zero seems to be consistent with SU(6) symme- 



try 



and recent lattice QCD results 27 



2 



should be sufficient to consider for it a dipole depen- 

-, where one would expect 



dence, C*{q 2 ) = {1 _ q y Ml ^ 
Maa ~ 0.85 — 1 GeV, to guarantee an axial transition ra- 
dius 2 R A in the range of 0.7 - 0.8 fm, and C^(O) - 1.2, 
which is the prediction of the off-diagonal Goldberger- 

Treiman relation (GTR), C^(0) = J^U^ = 1-2, with 
the irNA coupling /* = 2.2 fixed to the A width and 
/,r ~ 93 McV, the pion decay constant. 

There is no constraint from and lattice calcula- 

tions are still not conclusive about the size of possible vio- 
lations of the GTR. For instance, though values for (0) 
as low as 0.9 can be inferred in the chiral limit from the 
results of Ref. [2% they also predict C^(0)/ {yjlf^ 
to be greater than one. 



A second re- analysis 23j brings in the discussion two 
interesting points. First that both ANL and BNL data 
were measured in deuterium, and second, the uncertain- 
ties in the neutrino flux normalization. Deuteron struc- 
ture effects in the vd — > fi~A ++ n reaction, sometimes 
ignored, were estimated from the results of Ref. [l8[ to 
produce a reduction of the cross section from 5-10%. In 
what respects to the ANL and BNL flux uncertainties, 
the procedure followed in [23| is not robust from the sta- 
tistical point of view, since it ignores the correlations of 
these systematic errors 4 . Nevertheless, this latter work 
constitutes a clear step forward, and from a combined 
best fit to the ANL & BNL data, the authors of [HI 
find Cf{0) = 1.19 ± 0.08 in agreement with the GTR 
estimate. 



C£{q 2 ) ASSUMING AP DOMINANCE 

Traditionally, Adler's model and the GTR have been 
assumed, being the Maa axial mass adjusted in such 
a way that the AP contribution alone would lead to a 
reasonable description of the shape of the BNL q 2 dif- 
ferential v^p — > fi~pir + cross section (see e.g. Ref. 
These fits also describe reasonably well the q 2 depen- 
dence of the ANL data and the BNL total cross section 
but overestimate the size of the ANL data by 20% near 
the maximum [20]. Thus, ANL data might favor C^(0) 
values smaller than the GTR prediction. 

Recently, two re-analysis have been carried out trying 
to make compatible the GTR prediction for C^(0) and 
ANL data. In Ref. [Hj], Cg(0) is kept to its GTR value 
and three additional parameters, that control the C^(q 2 ) 
fall off, are fitted to the ANL data. In fact C^(q 2 ~ 0) 
is not so relevant due to phase space, and what is actu- 
ally important is the C${q 2 ) value in the region around 
—q 2 ~ 0.1 GeV 2 . Although ANL data are well repro- 
duced, we find the outcome in [22| to be unphysical, be- 
cause it provides a quite pronounced q 2 — dependence that 
gives rise to a too large axial transition radius 3 of around 
1.4 fm. Moreover, neither the fitted parameter statisti- 
cal errors, nor the corresponding correlation coefficients 
are calculated in (22[. Undoubtedly, the fit carried out 
there should be quite unstable, from the statistical point 
of view, because of the difficulty of determining three pa- 
rameters given the limited range of q 2 values covered in 
the ANL data set. Furthermore, the consistency of these 
results with the BNL data has not been tested. 



AXIAL FORM FACTORS INCLUDING THE 
CHIRAL NON-RESONANT BACKGROUND. 



2 It is denned from C£ (q 2 ) /Cf (0) = 1 + q 2 R 2 A /6 + 0(q 4 ). 

3 Further details and possible repercussions in neutrino induced co- 
herent pion production calculations are discussed in |2S(|. There, 
ANL data fits of the type proposed in [22j . but including chiral 
non-resonant contributions are also performed, finding that then 
the axial transition radius becomes even larger, about 2.5 fm. 



All the above-mentioned determinations of C§(q 2 ) 



suf- 



fer from a serious theoretical limitation. Though the AP 
mechanism dominates the neutrino pion production re- 
action, specially in the A ++ channel, there exist sizable 
non-resonant contributions of special relevance for low 
neutrino energies (below 1 GeV) of interest in T2K and 
MiniBooNE experiments. These background terms are 
totally fixed by the pattern of spontaneous chiral sym- 
metry breaking of QCD, and are given in terms of the 
nucleon and pion masses, the axial charge of the nucleon 
and the pion decay constant. When background terms 
are considered, the tension between ANL data and the 
GTR prediction for C^(0) substantially increases. In- 
deed, the fit carried out in [21] to the ANL data finds 
a value for C^(0) as low as 0.87 ± 0.08 with a reason- 
able axial transition radius of 0.75 ± 0.06 fm, and a large 
Gaussian correlation coefficient (r = 0.85), as expected 
from the above discussion of the results of Ref. [22j . 

Here, we follow the approach of Ref. [2l[, but imple- 
menting four major improvements: i) we include in the 
fit the BNL total v^p — > yTp-K + cross section measure- 
ments of Ref. 15| . Since there is no cut in the outgoing 



4 There exist some other aspects that might require further in- 
vestigation. For instance, additional parameters Panl an d Pbnl 
are introduced in |23| (see x 2 function in Eq. (37)) to account for 
the flux uncertainties. At very low q 2 values, da/dq 2 is totally 
dominated by C?g . If we had infinitely precise statistical mea- 
surements, the fit carried out in [23j would provide a very precise 
determination of the ratio C^-(0)/ \/p, but not of the form factor 
C^{0) . However, in such situation, one expects to extract (7^(0), 
though with an uncertainty dominated by that of the neutrino 
flux normalization. Besides, the fit to the BNL data uses the 
total cross-section data, for which the hadronic invariant mass is 
unconstrained, and the neutrino energy varies in the range 0.5—3 
GeV. Above 1 GeV, heavier resonances than the A(1232), and 
not considered in |23| . should play a role [24(1 . 
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pion-nucleon invariant mass in the BNL data, and in or- 
der to avoid heavier resonances from playing a significant 
role, we have just included the three lowest neutrino en- 
ergies: 0.65, 0.9 and 1.1 GcV. We do not use the BNL 
measurement of the q 2 — differential cross section, since 
it lacks an absolute normalization, ii) we take into ac- 
count deuteron effects in our theoretical calculation, Hi) 
we treat the uncertainties in the ANL and BNL neutrino 
flux normalizations as fully correlated systematic errors, 
improving thus the treatment adopted in Ref. 23j, and 
finally iv) in some fits, we relax the Adler's model con- 
straints, by setting C£ A {q 2 ) = C'£ A (0) (^(g 2 )/C 5 »), 
and explore the possibility of extracting some direct in- 
formation on (0). 

Let us consider first the neutrino-deuteron reaction 
vd — > fi~pTr + n measured in ANL and BNL. Owing to the 
inclusion of background terms, the formalism of Ref. (lij . 
where the pir + pair was replaced by a A ++ , cannot be 
used to account for deuteron corrections, and we must 
work with four particles in the final state. Neglecting 
the D— wave deuteron component and considering the 
neutron as a mere spectator, we find for the differential 
cross section on deuteron 



da 



dq 2 dW 



= J d 3 Pd \y d (p d )\ 2 



M da 



E p4 dq 2 dW 



p— off shell 



(i) 



where E p _ d — ra d — yf M 2 + p 2 , with m d the deuteron 
mass, is the energy of the off-shell proton inside the 
deuteron which has four-momentum p^ = (E p ^ d ,p d ). W 
is the final pir + invariant mass. The differential cross 

section da/dq 2 dW is computed using the model 

p— offshcll 

of Ref. [21] . Finally, ^ is the S— wave Paris potential 
deuteron wave function [29j normalized to 1. 

In what respects to the neutrino flux normalization un- 
certainties, we consider them as sources of 20% and 10% 
systematic errors for the ANL and BNL experiments re- 
spectively (see discussion in 23] ) . We have assumed that 
the ANL and BNL input data have independent statisti- 
cal errors (<7j) and fully-correlated systematic errors (ej), 
but no correlations linking the ANL and BNL sets. We 
end up with a 12 x 12 covariancc matrix, C, with two 
diagonal blocks. The first 9x9 block is for the ANL 
flux averaged q 2 — differential I'd ^p-n + n cross sec- 
tion data (with a 1.4 GeV cut in W), while the second 
3x3 block is for the BNL total cross sections mentioned 



above. Both blocks have the form C; 



The x 2 function is constructed by using the inverse of 
the covariance matrix. 

Results from several fits are compiled in Table U from 
where we draw several conclusions. First, by comparing 
fit II* with Ref. 2l|| . we deduce that the consideration 



of BNL data and flux uncertainties increases the value of 
C^(0) by about 9%, while strongly reduces the statisti- 
cal correlations between Cg (0) and Maa- Second, the 
inclusion of background terms reduces C5 (0) by about 



13%, while deuteron effects increase it by about 5%, con- 
sistently with the results of 21 1 and 3, 23[, respectively. 
Third, the fitted data are quite insensitive to 634(0), as 
fit V-VII results show. This is easily understood, taking 
for simplicity the massless lepton limit. In that case 



da 
dq^ 



^{[C£m 2 + q 2 a(q 2 )} 



(2) 



and 634(0) start contributing to a(q 2 ), i.e. to 0(q 2 ), 
which also gets contributions from vector form factors 
and terms proportional to dC^/dq 2 1 ?2 _ Q . This also ex- 
plains the large statistical correlations displayed in fits 
V-VII. Moreover, dC£ A /dq 2 \ 2 _ Q appears at order 0(q 4 ), 

which has prevented us to fitting the q 2 — shape of these 
form factors. Fourth, fit IV is probably the most robust 
from the statistical point of view. In Fig. [1] we display 
fit IV results for the ANL and BNL I'd —> ^Tp-K + n cross 
sections. Looking at the central values of C^(0), we con- 
clude that the violation of the off-diagonal GTR is about 
15% smaller than that suggested in Ref. 21], though it 
is definitely greater than that claimed in [23j , mostly be- 
cause in this latter work background terms were not con- 
sidered. However, GTR and fit IV C^(0) values differ in 
less than two sigmas, and the discrepancy is even smaller 
if Adler's constraints are removed. These new results 
are quite relevant for the neutrino induced coherent pion 
production process in nuclei which is much more forward 
peaked than the incoherent reaction. For instance, we 
expect the results in Ref. (30| . based in the determina- 
tion of C,f (0) of Ref. |21|, to underestimate cross sections 
by at least 30%. 

By using a state of the art theoretical model, we have 
determined the NA axial form factors from statistically 
improved fits to the combined ANL & BNL data. The in- 
clusion of chiral background terms significantly modifies 
the form factors. We have found violations of the GTR 
at the level of 2a, when the usual Adler's constraints 
are adopted. This will influence background calculations 
for T2K and MiniBooNE low energy neutrino precision 
oscillation experiments. 
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FIG. 1: Comparison of the ANL da/dq 2 differential (left panel) and ANL & BNL total (right panel) cross section data with fit IV 
theoretical results. Theoretical 68% confidence level bands are also displayed. Data in both plots include a systematic error (20% for ANL 
and 10% for BNL data) added in quadrature to the statistical ones. In the left panel, both data and results include a cut W < 1.4 GeV. 
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